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ABSTRACT 

Two dimensional quantum gravity coupled to a conformally invariant matter 
field of central charge c = n/2, is represented, in a discretized version, by n inde- 
pendent Ising spins per cell of the triangulations of a random surface. The matrix 
integral representation of this model leads to a diagrammatic expansion in powers 
of the cosmological constant for fixed genus. From the behaviour of this expansion 
at large orders, when the Ising coupling constant is tuned to criticality, one extracts 
the values of the string susceptibility exponent. We extend our previous calculation 
to order eight for genus zero and investigate now also the genus one case in order to 
check the possibility of having a well-defined double scaling limit even for c > 1. 
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A discretized approach to two-dimensional quantum gravity coupled to matter 
fields of central charge c = n/2, is easily formulated as an integral over 2 n N x N 
Hermitian matrices' 11 . In this model each cell of a discretized random surface 
carries n Ising spins. For a fixed triangulation these n species do not interact and 
the matter partition function is simply the n-th power of one Ising spin; however 
the various species are effectively coupled by the sum over the triangulations. In 

[2] 

a previous article , we found that at relatively small order in perturbation theory 
in the cosmological constant, we recovered with a surprizing accuracy, the known 
results for the string susceptibility exponent 7 string, for c < 1 and planar surfaces. 
For larger values of c, at least up to c = 4, we did not see any sign of the tachyonic 
instability of the theory. Similar conclusions on the absence of instabilities were 
obtained through Monte-Carlo simulations of the same model' 31 , but the values 
of the exponents were not reported. There were several attempts to evaluate the 
string susceptibility exponent for c > 1, by generating random surfaces embedded 
in a space of more than one dimension, but the results have not really converged 
so far' 41 . 

In our previous work , we had expanded the free energy up to sixth order in 
the cosmological constant g, for surfaces of genus zero. In this letter we extend 
the series to order eight, and include surfaces of genus one, in order to check the 
variation of the string susceptibility with genus which governs the existence of a 
double scaling limit. The expansion in powers of g is given by 4 Feynman diagrams 
in which each vertex carries n spins, i.e. 2 n "flavors". The Ising interaction (3 is 
present in the relative weight of the up-up and up-down propagators 161 through the 
parameter a defined as 



The sum over the n-Ising species is, for a given diagram, the n-th power of a 
one Ising sum, which is a polynomial in the variable a 2 . When the parameter a 
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varies and reaches its critical value, for instance 1/4 for one Ising spin, the theory 
crosses over from pure gravity, with '-/string — —1/2 to a coupling with conformally 
invariant matter with central charge n/2 and an exponent 7 string, which is for one 
Ising spin, but that we have tried to determine for n spins. 

The vertex of the 4 hermitian matrix model, carries four different lines of 
different colors i, % — 1,...,N. At Z-th order the permutations of the 4/ color 
indices, generate the various diagrams, and this requires a fast algorithm. If we 
use unitary instead of hermitian matrices, we obtain the same type of diagrams, 
with two incoming and two outgoing lines at each vertex, and the diagrams differ 
only by symmetry factors. In fact these symmetry factors are the same for the two 
ensembles at genus zero, but there are differences at genus one. With this unitary 
ensemble we generate the Feynman diagrams with the same tree sorting algorithm 
that was used for the study of the perturbation theory at large order of the Landau- 
Ginzburg Hamiltonian . In this method, the symmetry factor is automatically 
taken acount by the generation of all possible line combinations. Disconnected 
diagrams are avoided by a computation of the determinant of the adjacency matrix; 
it is made of 2 on the diagonal , and -1 for the elements corresponding to vertices 
connected by a line and some of the -1 of the element is replaced by 0. This 
determinant is equal to the number of Euler trails for a connected graph, but it 
vanishes for a disconnected diagram.' 81 

We have a certain number of checks on our eighth order result: the one Ising 
spin first' 91 . Then, if either a or n vanish we must recover the series of pure 
gravity' 101 . The free energy is given by 

^ = E^ fe ( 2 ) 

k 

and the coefficient behaves for the large order as 

C k ~ A k k^ +lst "^ (3) 
The ratio fa = Cfc/Cfc-i gives the value of the string susceptibility exponent by the 
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extrapolation, for example we have used the following ratio method in the previous 
analysis.' 111 

A k = (l + k)f k - kj k -i (4) 

1 string = 3 - fc(l + fc)(/fc - fk-l)/A k (5) 

Another method of the analysis, which we employed before, is the Pade approxi- 
mation for the ratio f k by 

1 (14- —) 

from which the string susceptibility ^string = 2 + b\ — ci is estimated. The unknown 
values A, b\ and c\ are determined from ftJk-i an d fk-2 f° r the fixed order k. 
This Pade becomes exact in the pure gravity a = case with the value b\ = —1/2 
and c\ = 2. This analysis is more convergent and it is convenient to observe the 
singularity due to the finite order perturbation, especially for the separated peaks 
at the large value of a. This method however fails to obtain the correct result for 
n = lie = 1) since (6) has no logk correction. For n = 2, we have checked that 
our eighth order perturbation agrees with — A k k~ s /logk behavior, for a = 1/4. 
So indeed we have 7 string = for n = 2, although we could not determine with 
precision the critical value 1/4 of a. The 1/logk in C k is indeed compatible with 

[121 

the logarithmic deviation to KPZ scaling found in the c=l problem. 

In Fig.la-lc, the string susceptibility ^string = 2 + b\ — c\ in (6) is plotted for 
n=l, n=4, and n=6 respectively. At a = and a — 1, the string susceptibility 
exponent becomes exactly 7 string = —1/2. And it takes the crossover from the pure 
gravity fixed point to the new fixed point at the critical vaue of a. In Fig.la-lc, 
the lines of a,b,c and d correspond to the Pade approximaion (6) for the string 
susceptibility exponent of order g 5 ,g®,g 7 and g 8 respectively. In Fig la (n=l), we 
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find that at the critical point a = 1/4, '-/string = —0.34, which is very close to the 
exact value of -1/3. The critical value of a is known from the exact solution.' 131 For 
n = 3, we have a slightly negative peak value of 7 string by this method. There 
is a logarithmic singularity for n=2, this may still influence the n=3 series, even 
if there is no log singularity in that case. The ratio metod of (5) for n = 3 case 
gives ^string — 0.02. For n — 4, ^string becomes 0.04 as shown in Fig lb. The n=2 
logarithmic singularity may still lead to underestimating the exponent for n=4. In 
the case n = 6 (Fig.lc), a second peak appears for larger values of a, but its height 
decreases at higher orders, and this second peak seems to merge with the first one 
near a = 0.2. We estimate the string susceptibility '"/string as 0.14 for n=6. The 
broken line is a guide line and the crossing point is the extrapolated value. In n=8, 
the second peak becomes more pronounced but it also approaches to the first small 
peak around a = 0.15. In this case, assuming the two peaks merge at the same 
point, we estimate ^string — 0.25. The obtained values up to n=8 (c=4) are not 
far from (c — 1)/12, which is the real part of the KPZ formula, 1141 valid for c < 1, 



_c-l y/(c-l)(c-25) 

Istring ~ yi 12 

although our values are lower. 

The free energy of the genus one is evaluated from the diagram of the genus 
zero by permutation of the four lines at the vertices. We calculate the ratio of 
each order of genus for a fixed diagram, and then normalize the coefficients of each 
diagram by the genus zero results. The genus one case for the pure gravity has 

[15] 

been studied , from which we checked our perturbation series. The large order 
behavior of the coefficient Ck of the genus one of the free energy is also written by 



C k ~ A k k- 3+ ^s (8) 
We have a scaling relation of ^string for /i-genus case, 
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Istring — Istring + ^Istring (9) 



For the genus one, we have h — 1, 



Istring + Istring ~ ^ (^) 



if c < 1 [161 . This relation may be violated for c > 1 and we are unable to assume 
this relation. The free energy of the genus one for the Hermitian matrix model is 
expanded up to order six, and here we present the result up to fifth order result, 
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Ei(g) = -g + <? 2 [20(1 + a 2 ) n + 10(1 + a 4 )'"] 

- /[256(1 + a 2 ) 2n + 416(1 + a 2 + 2a 4 ) n + ^(1 + 3a 2 ) n + ^(1 + 3a 4 ) n ] 
+ /[2816(1 + a 2 f n + 6400(1 + 4a 2 + 3a 4 ) n + 2080(1 + 6a 2 + a 4 ) n 

+ 896(1 + a 2 + 4a 4 + a 6 + a 8 ) n + 880(1 + 6a 4 + a 8 ) n 

+ 9856(1 + a 2 + 5a 4 + a 6 ) n + 6144(1 + 2a 2 + 5a 4 ) n + 4480(1 + 3a 2 + 3a 4 + a 6 ) n 
+ 3712(1 + 2a 2 + 3a 4 + 2a 6 ) n + 2912(1 + 5a 4 + 2a 6 ) n ] 

- /[28672(1 + a 2 ) 4n + 63488((1 + a 2 )(l + 4a 2 + 3a 4 )) n 

+ 32768((1 + a 2 )(l + 2a 2 + 3a 4 + 2a 6 )) n + 78848((1 + a 2 )(l + 6a 2 + a 4 )) n 
+ 77824((1 + a 2 )(l + 3a 2 + 3a 4 + a 6 )) n 
+ 107520((1 + a 2 )(l + 2a 2 + 5a 4 )) n 

+ 81920((1 + a 2 )(l + a 2 + 5a 4 + a 6 )) n + 34816(1 + 3a 2 ) 2n 

+ 34304((1 + 3a 2 )(l + a 2 + 2a 4 )) n + 21606.4(1 + 10a 2 + 5a 4 ) n 

+ 49664(1 + 6a 2 + 5a 4 + 4a 6 ) n + 125952(1 + 4a 2 + 9a 4 + 2a 6 ) n 

+ 92672(1 + 3a 2 + 7a 4 + 5a 6 ) n + 47616(1 + 2a 2 + 7a 4 + 4a 6 + 2a 8 ) n 

+ 18944(1 + 3a 2 + 6a 4 + 3a 6 + 3a 8 ) n + 176128(1 + 2a 2 + 9a 4 + 4a 6 ) n 

+ 60928(1 + a 2 + 9a 4 + 3a 6 + 2a 8 ) n 

+ 38912(1 + 3a 2 + 11a 4 + a 6 ) n + 97792(1 + a 2 + 7a 4 + 7a 6 ) n 
+ 44032(1 + 2a 2 + 10a 4 + 2a 6 + a 8 ) n + 102400(1 + a 2 + 8a 4 + 5a 6 + a 8 ) n 
+ 6400(1 + 2a 2 + 5a 4 + 4a 6 + 4a 8 ) n + 29696(1 + a 2 + 7a 4 + 3a 6 + 4a 8 ) n 
+ 4710.4(1 + 10a 4 + 5a 8 ) n + 29184(1 + 8a 4 + 4a 6 + 3a 8 ) n 
+ 46592(1 + 7a 4 + 6a 6 + 2a 8 ) n + 63488((1 + a 2 )(l + 4a 2 + 3a 4 )) n 
+ 3942.4(1 + 5a 4 + 10a 6 ) n ] + . . . 



[17] 

At a = 0, we have the correct expression obtained before, 



Ei(g) = -g + 30g 2 - 1056# 3 + 40176/ - 1600819.2# 5 
+ 65774592/ - 2762461769/ + 1.17944875 x 10 11 / 



(11) 



(12) 
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The string susceptibility ^string is obtained by the ratio method of (5) as shown 
in Fig. 2a and Fig. 2b for n—1 and n=3 respectively. Since the successive lines in 
Fig. 2 are approaching smoothly to the line ^string— ^ especially for the small value 
of a (a < 0.25), the results support the result ^string = 2. For n = 3 there is a tiny 
peak around a = 0.2. Even if this relation is violated, the value of ^/string remains 
close to two. Other cases of different n have similar curve as Fig. 2 and the relation 
of (10) seem to hold at least when the critical value a remains to be small enough. 

For n = 2 case, we find that the coefficient Ck in genus one behaves as (8) 
without log k correction. With the assumption of the log k correction, we are 
unable to hold the relation ^string = 2. This rules out logarithmic terms at genus 
one. 

We have also developed the expansion for the genus- two case up to order g 6 . 
By the simple ratio method, in which the ratio Ck/Ck-i is plotted against 1/k, 
we find the same value of the critical cosmological constant A from three different 
genus cases. It seems that genus one case gives most reliable value of the critical 
cosmological constant A. The scaling relation of (9) seems to hold. The renor- 

[18] 

malization group approach for the matrix model has been investigated. We will 
discuss this approach for the present model elsewhere. 

The authors thank Al.Zamolodchikov for his interest and his comments. 
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Figure caption 

Fig.la ^string for n=l Ising case. At a=l/4, ^string becomes -0.34. The lines a,b,c 
and d are Pade analysis of (6) from order g 5 ,g 6 ,g 7 and g 8 , respectively. 

Fig.lb '"/string for n=4 (c=2). a,b,c and d are same as Fig.la. 

Fig.lc ^string for n=6 (c=3). a,b,c and d are same as Fig.la. The broken lines are 
guided lines for the extrapolation. 

Fig.la. ^string of n=l for genus one. The broken line is ^string — 2. The lines 
a,b,c and d are the result of the ratio method (5) for order g 3 ,g 4 ,g 5 and g 6 , 
respectively. 

Fig.lb ^string of n=3 for genus one. The lines a,b,c and d are same as Fig.2a. 
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